Tuning the interlayer spacing in high Tc superconductors: penetration depth and 

two-dimensional superfluid density 
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Substantial control of the interlayer spacing in Bi-based high temperature superconductors has 
been achieved through the intercalation of guest molecules between the superconducting layers. 
Measurements using implanted muons reveal that the penetration depth increases with increasing 
layer separation while Tc does not vary appreciably, demonstrating that the bulk superfluid density 
is not the determining factor controlling Tc. Our results strongly suggest that the superfluid density 
appearing in the Uemura scaling relation n^jm* ex Tc should be interpreted as the two dimensional 
density within the superconducting layers, which we find to be constant for each class of system 
investigated. 
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A fundamental property of high temperature supercon- 
ductors (HTSs) is the presence of Cu02 planes with some 
degree of weak coupling between them [IJ. The nature 
of this interlayer coupling ^] has been fiercely debated 
in the past and recent advances in elucidating the Fermi 
surface topology of the cuprates [3|, |j, |5| make an under- 
standing of the effective dimensionality of the the HTSs 
a priority. In order to address the question of dimension- 
ality experimentally we have modulated the interlayer 
coupling of the Cu02 planes in two HTS systems by in- 
tercalating molecular groups between the layers. This al- 
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FIG. 1: (a) Schematic of the Bi2212 intercalates with molecu- 
lar groups acting to separated the Cu02 bilayers. (b) Exam- 
ple TF /i"'"SR spectra measured in 40 mT above and below Tc 
for polycrystalline Bi2212. The strong dephasing for T < Tc 
is caused by distribution of fields in the vortex lattice, (c) 
Interlayer/interbilayer separation d for each intercalate of the 
Bi2212 and Bi2201 series. 



lows us to approach the two dimensional (2D) limit with- 
out perturbing the superconducting oxide block. Here 
we present the results of a muon spin rotation (/x+SR) 
investigation of the evolution of the penetration depth 
(and hence ng/m*) with layer separation. We are able to 
demonstrate that the bulk superfluid density rig is not the 
parameter controlling Tc and that the superfluid density 
appearing in the Uemura relation ng/m* oc Tc should be 
interpreted as the superfluid density within the 2D su- 
perconducting layers, which is independent of the layer 
separation. 

The bismuth based cuprates possess weakly bound 
Bi2 02 double layers that allow the intercalation of guest 
molecules without any substantial change in the super- 
conducting block [i, 0, d, H, [l^, [iH. Weak interaction 
between the guest molecule and host layers ensures that 
the intercalation process does not distort the supercon- 
ducting layers. This overcomes a substantial difficulty 
with alternative attempts to control the interlayer spac- 
ing I4I , which involve the sequential deposition of super- 
conducting and insulating layers. Moreover, the interca- 
lation of long chain organic molecules (shown schemati- 
cally in Fig. [IJa)) [71, Isj allows the possibility of tuning 
the layer (or bilayer) separation d on a scale fixed by the 
length of the hydrocarbon chain. Organic guest species 
electrically isolate the layers since hydrocarbons are, in 
general, very good insulators. In this paper we exam- 
ine intercalation series based on two classes of Bi-based 
superconductors, namely bilayer Bi2Sri.5Cai.5Cu208-i-5 
(Bi2212) and single layer Bi2Sri.6Lan 4Cu06-t. (Bi2201). 
Previous work on these classes [6|, [Zl, |8[ Isl, llOl |ll| has 
suggested that while the intercalation of guest species al- 
lows d to be varied by a factor of three, the transition 
temperature of the intercalated materials does not vary 
appreciably compared to the pristine compound. 

Transverse field muon spin rotation (TF /i+SR) pro- 
vides a means of accurately measuring the internal mag- 
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FIG. 2: (a) Temperature dependence of the magnetic field 
distribution width B^ms for Bi2212 (left) and Bi2201 (right) 
intercalates with Ba = 40 mT. Fits are described in the main 
text, (b) The evolution of Brms with varying applied field at 
T = 5 K (Bi2201) and T = 20 K (Bi2212) shows the general 
trend of a rounded peak in applied field. Additional data for 
pristine Bi2212 were obtained from Ref . [ig. The results of the 
fits (described in the main text) are shown with a solid line 
for the range over which the fitting was carried out; extrapo- 
lations are shown with a dashed line. The materials with the 
largest layer separations show markedly different behavior in 
both temperature and applied field (dotted lines are guides 
to the eye). 



netic field distribution caused by the vortex lattice (VL) 
in a type II superconductor [13| . The vortex state of the 
pristine polycrystaUine Bi2212 and Bi2201 systems have 

Mm 



been studied with /x+SR previously [M, ll5|, 116| ■ In a TF 
/i^SR experiment spin polarized muons are implanted 
in the bulk of a superconductor, in the presence of a 
magnetic field Bd < B < Bc2, which is applied perpen- 
dicular to the initial muon spin direction. The muons 
stop at random positions on the length scale of the VL 
where they precess about the total local magnetic field 
B at the muon site (mainly due to the VL), with fre- 
quency LUfj_ = ^i_iB, where 7^ is the muon gyromagnetic 



ratio (= 27r x 135.5 MHz T~^). The observed property 
of the experiment is the time evolution of the muon spin 
polarization P^ (t) , which allows the determination of the 
distribution of local magnetic fields across the sample 
volume p{B), via 



/•OO 

P,(i) = / p{B) cos(7^Bi + <j>)dB, 
Jo 



(1) 



where </> is a phase offset associated with a particular 
detector geometry. 

Powder samples of several intercalates of Bi2212 and 
Bi2201 were prepared as reported previously 0,0,11,0]. 
The intercalates studied and the separation d of CUO2 
layers (or bilayers) achieved are given in Fig. [ijc) and 
Table m /x+SR measurements were made at the ISIS Fa- 
cility using the MuSR spectrometer and the S/iS facility 
using the GPS spectrometer. Samples were pressed into 
pellets and mounted on a hematite backing plate to re- 
duce the background signal. 

Fig. [ijb) shows example TF /i+SR spectra measured 
above and below Tc for pristine Bi2212 in an applied 
field of 40 mT. Above Tc some slight broadening of the 
spectrum is attributable to dephasing of the muon spins 
caused by randomly directed nuclear moments near the 
muon stopping sites. Below Tc the spectrum broadens 
considerably due to the dephasing contribution from the 
VL. The expected width of the field distribution in the 
powder sample was estimated using the numerical results 
reported for the field width in the Ginzburg-Landau (GL) 
model that were obtained for a field normal to the lay- 
ers ITJ. A polycrystaUine average was taken under the 
assumption that the length scales A and ^ diverge follow- 
ing l/cos0 as the field orientation approaches the plane 
at 6* = (high anisotropy limit), with the corresponding 
contribution to the overall width scaling as cos 6'. Fit- 
ting our data with these results allowed us to extract the 
second moment of the distribution, i^rms, at each mea- 
sured field and temperature, and relate it to the vortex 
properties and penetration depth of the material. 

Fig. ^a.) shows the temperature dependence of Brms 
in an applied field of 5 mT for each intercalate studied. 
For both series, the magnitude of i?rms falls continuously 
with increasing temperature to a low constant value for 
T > Tc- The exceptions are the materials with the largest 
layer separation d in each series (discussed below). Fits 
are shown to an empirical power law functional form [l6[ 
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where i?vL(0) is the zero temperature contribution from 
the VL and i?bg represents a background contribution 
from the dipolc fields from nearby nuclei. The parame- 
ters derived from these fits are given in table IH There 
is little variation in Tc while i?vL(0) decreases smoothly 
with increasing layer spacing. The curvature parameter 
r falls within the range 2 < r < 6, which is a slightly 



TABLE I: Properties of the intercalates of Bi2Sri,5Cai.5Cu208+i (Bi2212) and Bi2Sri.6Lao.4Cu06+£ (Bi2201), with parameters 
derived from the fitting routines described in the main text (Py=pyridine, Me=methyl) . Fitted values are given for Eq. ([2]) 
and the fitting routine described in the main text. Fits were not possible for the materials with largest d from each class, 
although the values of Tc for these materials (derived from magnetization measurements) were found to coincide with those of 
the respective series to within ~ 1 K0- 



Intercalant 


d (nm) 


T, (K) 


r 


Bvl(O) (mT) 


Bbg (mT) 


Aab (nm) 


C(nm) 


d/AL(xl0*m-i) 


Bi2212 


none 

[Hgl2]0.5 

[(Py-CH3)2Hgl4]o.35 

[(Me3S)2Hgl4]o.34 

[(Py-C8Hl7)2Hgl4]o.35 


1.53 
2.25 
2.61 
2.79 
3.82 


82.4(2) 

80.6(1) 

79.0(1) 

88(3) 


4.6(6) 
5.7(4) 
2.6(1) 
2.1(3) 


1.52(5) 
1.30(3) 
1.04(1) 
0.90(6) 


0.07(2) 
0.09(1) 
0.11(1) 
0.07(1) 


200(6) 
230(10) 

271(7) 
270(8) 


6(1) 
3(1) 




3.8(2) 
4.3(4) 
3.6(2) 
3.8(2) 


Bi2201 


none 

(Hgl2)0.5 
[(Py-CH3)2Hgl4]o.35 


1.21 
1.92 
2.30 


22.1(2) 
22.3(5) 


3.0(2) 
2.1(3) 


0.98(3) 
0.63(5) 


0.006(3) 
0.05(3) 


290(8) 
370(10) 


13(3) 
10(2) 




1.44(8) 
1.40(8) 



larger range than that previously observed when this ap- 
proach was applied to YBasCuaO^, (YBCO) 0. The 
background contribution Bw is within the range found 
for pristine materials IJ, |l6| , confirming that the inter- 
calant is not introducing significant broadening. 

The magnetic field dependence of Srms [l8| measured 
at low temperature for each series is shown in Fig.[2lb). 
The general trend is for i?rms to increase sharply with in- 
creasing field before describing a broad peak and decreas- 
ing for high fields. Again the exceptions to the general 
behavior occur for the largest layer separation d. 

To extract the penetration depth Aab from our mea- 
surements we compare the field dependence of i?rms with 
the behaviour expected for the GL model (Fig.[2ljb)) [l7|. 
The model describes the data well for fields around the 
peak value of i?rms and above, but provides a poor de- 
scription of the data at low fields (probably due to the 
anisotropic magnetization and the distribution of demag- 
netizing factors in the differently aligned crystallites). 
We note that the ^ dependence of the theoretical line 
width becomes very weak near the peak, making an accu- 
rate assignment of ^ unimportant in our discussion of Aab- 
Indeed, close to the expected peak in Srms, for k > 70 
(as is the case for our systems), the GL model predicts 



i3i, im. 



the well known result B^^,, = O.OOSTl^gA"'' 
(Here, the penetration depth A should be interpreted as 
the effective penetration depth which includes contribu- 
tions from supercurrents both parallel and perpendicular 
to the layers.) The values of Aab and ^ derived from our 
data are given in Table II|ig |. As above, we were not able 
to model the samples with the largest d for each class of 
material. 

The general trend revealed by our measurements is a 
clear increase in Aab with increasing layer separation d. 
In the clean limit, we expect A~^ = /ioe^ns/m*, where Ug 
is the bulk superfluid density; thus ng/m* decreases with 
increasing d, but crucially Tc remains almost unchanged. 
This strongly suggests that the bulk superfluid density 



is not the determining factor controlling Tc in these sys- 
tems. We note that our results differ from those involv- 
ing varying the thickness of a PrBa2Cu307 (PBCO) layer 
in PBCO/YBCO superlattices, where it was found that 
Tc oc 1/d [ij, [20|. We stress, however, that chemical 
intercalation should be expected to produce better iso- 
lated layers than those found in superlattice systems (see 
above). 

The behaviour we observe contrasts with that expected 
from a simple application of the Uemura scaling relation 
Tc oc Us/m* [2l| which holds for many underdoped super- 
conductors (Fig[3I^a)). This scaling relation is, however, 
often interpreted as suggesting the strong two dimension- 
ality of the superfluid in underdoped cuprate supercon- 
ductors [2^. We might, therefore, expect that for the in- 
tercalated systems Tc oc n^^ /m* , where nf^ is the super- 
fluid density in the superconducting layers. If we assume 
that the superfluid density in each Cu02 layer is the same 
for each class of material (i.e. a constant for Bi2212 and 
a different constant for Bi2201), then we should expect 
that d/Xl^ oc nf^/m* oc Tc. We see from Table U and 
Fig. ^h) that our data are quite well described by this 
assumption: d/X^^ is approximately constant for each 
class of materials, as is Tc- This strongly suggests that 
the superfluid density appearing in the Uemura relation 
should be interpreted as the 2D density in the supercon- 
ducting planes. 

Some insight into the meaning of our results may be 
obtained if we assume that underdoping destroys super- 
conductivity at a quantum critical point (QCP) in these 



systems. Near a QCP (23|, |2J| an energy scale such as 
Tc should vary as Tc oc 6^'^, where 5 is the difference in 
doping from the QCP, z is the dynamic exponent and 
h' the the correlation length exponent. Combining this 
with the prediction of Josephson scaling near a QCP, 
ns(0) oc S^'-+^-^>, we have Tc oc 71^(0)^/(^+^-2)^ p^j. 
D = 2 this predicts Uemura scaling Tc oc nf^{0) ^ . 
Our findings may therefore be suggestive of the influence 




FIG. 3: (a) Uemura plot of Tc versus 1/Aj,b for our com- 
pounds. The solid line represents the scaling proposed for the 
underdoped cuprates ^1|. Dotted lines are guides to the eye 
showing the observed departure from the Uemura scaling, (b) 
The variation of the penetration depth Aab with the inverse 
layer separation 1/d gives a straight line dependence for each 
series, suggesting the constancy of d/A^b(oc n^ /m*). 



of 2D quantum fluctuations in the underdoped regime. 

Finally we note that VL breakup effects explain the 
anomalous behavior of the samples with the largest layer 
separation (see Fig. 2). A sharp drop in Brnis{T), sim- 
ilar to that observed in our Bi2212 {d = 3.82 nm) and 
Bi2201 {d — 2.3 nm) samples, has previously been ob- 
served in pristine Bi2212 27|, where it was attributed to 
VL melting at small superfluid densities close to Tc- Both 
electromagnetic and Josephson coupling are expected to 



27|. The tem- 



stabilize the VL in these materials [26 

perature at which thermal fluctuations break up a VL 

stabilized solely by electromagnetic coupling is given by 



Tem = ^qc/ {32kBfJ.oT^^ ^ab) [28| (where c is the average 
Cu02 layer separation), which is constant for each se- 



ries as r ^ 0. We can therefore estimate Tc, 



14 K 



for Bi2212 and Tom -^ 10 K for Bi2201, both in reason- 
able agreement with the temperatures at which we ob- 
serve a sharp drop in Br^s{T) in the samples with the 
largest layer spacing. This suggests that the Josephson 
coupling, which decreases strongly with increasing layer 
spacing, stabilizes the VL to higher temperature in the 
samples with smaller layer spacing but is too weak to do 
so in the two samples with the largest d. This leaves only 
electromagnetic coupling and the observed VL breakup 
for T > Ten,. 
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